#5.1

esaligna = read.table("http://ecologia.ib.usp.br/bie5782/lib/exe/fetch.php?media=dados:esaligna.csv", sep = ",", header = T)
esaligna
grafico1 = plot(esaligna$dap~esaligna$ht, xlab = "", ylab = "", las=1, xlim=c(5,30), ylim=c(2,30), tcl= 0.8 ,bty="l",main="Grfico DAP x Altura")

mtext("Altura(m)", side=1, line=3, at=15,cex=1, family="mono") 
mtext("DAP(cm)", side=2, line=2.5, at=15,cex=1, family="mono") 


#5.2

par(mfrow = c(1,2))
dap = esaligna$dap
talhao = esaligna$talhao
box=boxplot(dap~talhao, data=esaligna ,xlab="talho", ylab="dap")

dap.m = aggregate(x=esaligna$dap,by=list(esaligna$talhao),FUN=mean) 
dap.sd = aggregate(x=esaligna$dap,by=list(esaligna$talhao),FUN=sd) 

unique(esaligna$talhao)
eixopos = dap.m$x + dap.sd$x
eixoneg = dap.m$x - dap.sd$x

plot(dap.m, las=1, xlab = "talho",ylab="dap", main = "media dap para cada talhao", ylim= c(0,30))
segments(x0=dap.m$Group.1, y0 = eixopos, y1 = eixoneg)

table(esaligna$talhao)


matlines(dap.m, lty = c(3,2), col = "red") 
mtext("A", side=2, line=42, las=1, at=33,cex=3, family="mono", col = "red") 
mtext("B", side=2, line=2,las=1, at=33,cex=3, family="mono", col = "red") 
x11

#5.3

par(mfrow = c(1,2))
adv = read.table("http://ecologia.ib.usp.br/bie5782/lib/exe/fetch.php?media=bie5782:01_curso2009:material:exercicio3.csv", sep = ",", header = T)

plot(adv$x1, adv$y1, pch=17, ylab = "Euclidean distances", xlab= "Log(Patch size)(ha)", ylim= c(0.0, 3.0), tcl=0.2)
segments(x0=0.5, y0=2.25, x1=2.5,  y1=1, lwd=2)
text(2.5,3.0,"a",cex=2,family="serif")

par(tcl=0.2)
boxplot(adv$x2~adv$y2, xlab="", ylab="", ylim = c(0.0, 3.0), xaxt="n", outline = F)
text(6.5, 3.0,"b",cex=2,family="serif")

axis(side = 1, labels = F)
mtext("Small", side=1, line=1, las=1, at=1,cex=1) 
mtext("Medium", side=1, line=1, las=1, at=2,cex=1) 
mtext("Edge", side=1, line=2, las=1, at=2,cex=1) 
mtext("Medium", side=1, line=1, las=1, at=3,cex=1) 
mtext("Interior", side=1, line=2, las=1, at=3,cex=1)  
mtext("Large", side=1, line=1, las=1, at=4,cex=1) 
mtext("Edge", side=1, line=2, las=1, at=4,cex=1)
mtext("Large", side=1, line=1, las=1, at=5,cex=1) 
mtext("Interior", side=1, line=2, las=1, at=5,cex=1) 
mtext("Control", side=1, line=1, las=1, at=6,cex=1)

text(1, 3.0,"*",cex=1,family="serif")
text(2, 3.0,"*",cex=1,family="serif")
text(3, 3.0,"**",cex=1,family="serif")
text(4, 3.0,"*",cex=1,family="serif")
text(5, 3.0,"***",cex=1,family="serif")
